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, | Abstract. Two results about equidistribution of tile orientations in primitive substitution 

tilings are stated, one for finitely many, one for infinitely many orientations. Furthermore, 
consequences for the associated diffraction spectra and the dynamical systems are discussed. 

1. Substitution tilings 

Several important models of quasicrystals are provided by (decorations of) nonperiodic 
substitution tilings. In this paper, only tilings of the plane are considered. In larger dimensions 
some aspects become rather complicated, in particular, generalisations of Theorem 12.21 A 
i -^h ■ substitution (or tile-substitution, but we stick to the shorter term here) is just a rule how to 

enlarge and dissect a collection of tiles Pi, ... , P m into copies of Pi, ... , P m , as indicated in 
Figure [1] (left) or Figure [2j Iterating this rule yields tilings of larger and larger portions of 
the plane. However, the precise description of a substitution tiling in M? is usually given as 
follows. 

Let Pi, ... , P m be compact subsets of M 2 , such that each Pj equals the closure of its interior. 
These are the prototiles, the pieces our tiling is built of. We refer to sets which are congruent 
to some prototile simply as tiles. (Sometimes two prototiles are allowed to be congruent, 
■ compare Figure [TJ Then we equip each one of them with an additional attribute (colour, 

, decoration) to distinct them. Whenever we speak of 'congruent' tiles in the sequel, this 

means 'congruent with respect to possible decorations'.) A tiling (with prototiles Pi, ... , P m ) 
is a collection of tiles (each one congruent to some Pj) which covers M? such that every i£l 2 
is contained in the interior of at most one tile. Note, that each such tiling T can be described 
by placements of the prototiles: 

^! T = {PaiPii + ti, R m Pi 2 + *2, • • •}, 

where R ai denotes a rotation through a^, and fj G K 2 are translation vectors. For later 
purposes, we define the orientation 4>(T) of a tile: If T = R a Pi+t, a £ [0, 2ir[, then (f>(T) = a. 

Now, let A > 1 and a be a map which maps each prototile Pi to a non-overlapping collection 
of tiles, each one congruent to some prototile Pj. Furthermore, we require the substitution 
a to be selfsimilar in the sequel. That is, APj = Ureo-fPi) ^ or ^ — * — m - By setting 
a(RPi + 1) = Rcr(Pi) + At, the substitution a extends naturally to any tile congruent to some 
prototile. In a similar way, by 0"({Ti, T%, . . .}) = <7(Xi)U<r(T2)U . . ., the substitution a extends 
to all (finite or infinite) collections of tiles. In particular, one can iterate a on the prototiles 
to obtain the A:-th order supertiles a k (Pi). 

Definition 1.1. Let a be a tile- substitution with prototiles Pi, . . . ,P m . A tiling T is called 
substitution tiling (with substitution a) if for each finite subset T C T there are i, k such that 
a copy of T is contained in some supertile a k (Ti). The family of all substitution tilings with 
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substitution a is denoted by X CT . The matrix M a with (M CT )jj = #{T 6 o~(Pj) '■ T = RPi + t} 
is called substitution matrix. A substitution is primitive, if there is k > 1 such that is 
strictly positive. 

Here, j^A denotes the cardinality of the set A. This definition of a substitution tiling fits 
well into several contexts, for instance, to the hull of a tiling, see Section [H 

2. Statistical circular symmetry 

Most of the well-known substitution tilings have the property that all prototiles occur in 
only finitely many orientations. For instance, the Penrose dart and kite tiling has just two 
prototiles, a convex and a non-convex quadrangle. Each one occurs in one of 10 different 
orientations throughout the tiling. There are also substitution tilings where the tiles occur in 
infinitely many orientations. The most prominent examples are certainly the pinwheel tilings 
|14j . But there are many other examples, compare [6], [7]. Two examples are shown in Figure 

El 

It was known for the pinwheel tiling that there are not only infinitely many orientations of 
the tiles, but that they are furthermore equidistributed [2], [12]. Intuitively, this means that 
each orientation occurs with the same frequency throughout the tiling. In order to make this 
precise, we need the following definition. Recall that a sequence (aj)j>o is equidistributed in 
[0, 2tt[, if for all < x < y < 2vr holds: 
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Since the above sum is not absolutely convergent, it depends strongly on the ordering of the 
sequence. Thus we need to order the tiles in the tiling in some not-too-weird way. This is 
done in the following definition. 

Definition 2.1. Let T = \T±,T2, ■ ■ ■ } be a primitive substitution tiling. Let the numbering 
be such that the sequence (Tj)j >1 satisfies the following property: for all n> I, there is some 
£ > n such that the patch {T\ . . . , TJ?} is congruent to some supertile a k {Pi) for some prototile 
Pi, k > 1. The tiling T a has statistical circular symmetry, if, for all < x < y < 2tt, one 
has: 



lim - V 1[ y] (0(Tj)) = 



The ordering in the definition is compatible with other natural notions, for instance, order- 
ing the tiles with respect to their distance to the origin. The following result was obtained in 
[6]. Only while preparing the present text, the author realized that a similar result already 
appeared in p3]; see also the comments in Section SJ 

Theorem 2.2. Let T be a primitive substitution tiling in M?, where copies of some prototile 
occur in infinitely many orientations. Then T is of statistical circular symmetry. 

In the next section this theorem is used to show that any such tiling has circular diffraction 
spectrum. Now we show an analogue of the theorem above for tilings with tiles in finitely 
many orientations. One can ask whether there are substitution tilings where tiles occur more 
frequent in one direction than the other. This is possible, see Figure Q] for an example. But 
this is achieved by substituting the horizontal rectangle differently from the vertical rectangle. 
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Figure 1 . A substitution tiling where horizontal rectangles are more frequent 
in the tiling than vertical ones. By Theorem 12.31 this can be achieved only 
when the two types of rectangles are substituted differently from each other. 



In fact we have two prototiles in this example, each one occurring in one direction only. This 
is illustrated in the right part of Figure [1] by assigning different decorations to the tiles. The 
following theorem states that each prototile Pi occurs in each of its orientations with the 
same frequency in a primitive substitution tiling. In order to make this precise, let freq(Pj, a) 
denote the frequency of tiles R a Pi + 1 in a tiling T. That is, 

m f , p ^ r #{T<=TnB r : 3t T = R Q P i + t} 
(1) freq(P,a,T)= hm #{T e T n B r} ' 

where B r denotes the closed ball of radius r around the origin. For primitive substitution 
tilings, this frequency exists and is independent of the choice of the origin. This can be seen 
as follows: Up to now, we identified types of tiles with respect to congruence congruence 
classes. In a tiling with finitely many orientations, one can as well identify types of tiles 
with respect to translation classes. Usually, this results in a larger number of prototiles (e.g., 
the Penrose dart and kite tilings have two prototiles up to congruence, but 20 prototiles 
up to translation). The frequency of each prototile in a primitive substitution tiling can 
be obtained from the substitution matrix M ff . This result is folklore, see for instance [13] , 
[5]: The normed Perron- Frobenius eigenvector of M a contains the relative frequencies of the 
prototiles. In particular, in all tilings in X CT , these frequencies are the same. Therefore, for 
any given primitive substitution, we can drop the T in the definition of the frequency in this 
case and write just freq(Pj, a). 

Theorem 2.3. Let a be a primitive substitution with prototiles {Pi, ■ ■ ■ , P m }- If freq(Pj, a) ^ 
0, freq(P, (3) / then freq(P;, a) = freq(P, (3) . 

Proof. If freq(Pj, a) ^ 0, then the prototiles occur in finitely many orientations only, by prim- 
itivity of the substitution and Theorem [221 Therefore we obtain finitely many prototiles with 
respect to translations: (Pi, 0), (Pi, a;i,i), (Pi, 01,2), • • • , (Pi, «i,fc); (P2, 0), (P2 

Consider a tiling T € X CT . By the remarks preceding the theorem, freq(Pj,0) is well 
defined, and its value is the same in the tiling T and in the tiling R ai T G X CT . And trivially, 
freq(Pj,0) in T equals freq(Pj, ctij) in R aij T. Therefore, freq(Pj,0) = freq(Pj, ajj), and the 
claim follows. □ 
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Figure 2. Two substitutions yielding tiling with statistical circular symme- 
try. The left one is due to Danzer and Goodman-Strauss [8], the right one is 
from [6]. 



3. Diffraction 

There is a wealth of literature devoted to the mathematical description of the diffraction 
spectrum of nonperiodic structures, see for instance [1], [9], [10], [18], or the introductory text 
[lj and references therein. Usually one starts with a discrete point set A C K 2 rather than 
a tiling. In the present context of tilings, this can be achieved by decorating each prototile 
with one or more points ('atomic decorations'). 

Let A be such a point set, obtained from a tiling. The diffraction spectrum of A is the 
Fourier transform of the autocorrelation of A. The autocorrelation (also called Patterson 
function) of A is 

( 2 ) 7 = 7A = lim — V 5 x - y , 

r^oo vol (B r ) 

if this limit exists. Here and in the following we are working in the framework of tempered 
distributions [17], [19]. Thus the limit above is to be understood as a vague limit. If the limit 
does not exist, there is a subsequence converging to an autocorrelation which is known to be 
translation bounded (since we started with a translation bounded measure 5\) and positive 
definite. Thus 7 can be regarded as a translation bounded positive definite measure. Now, 
the Fourier transform 7 of 7 is the diffraction measure, and its support is the diffraction 
spectrum of A. We are interested in the nature of 7. By Lebesgue's decomposition theorem, 
7 has a unique decomposition 

7 = %p + lsc + lac 

with respect to the Lebesgue measure in R 2 . The pure point part j pp is a countable sum 
of Dirac measures. The absolutely continuous part 7 ac is a measure with locally integrable 
density function supported on a set of positive Lebesgue measure. The singular continuous 
7 SC part assigns zero to all finite sets, but is supported on a set with Lebesgue measure zero. 
The (idealized) diffraction spectrum of a quasicrystal is pure point, the singular parts vanish 
completely. 
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The following results show that substitution tilings with statistical circular symmetry have 
a diffraction measure which is of (perfect) circular symmetry. Consequently, the pure point 
part of the diffraction is trivial. 

Theorem 3.1. Let T = {R^P^ + t\,R a , 2 Pi 2 + t2, ■ ■ ■ } be a primitive substitution tiling of 
statistical circular symmetry, with prototiles P%, . . . , P m . Choose 'control points' Xi £ Pi, 
and define A = {R ai Xi x + t\, R a2 Xi 2 + t%, . . .}. Then the autocorrelation 7a of A is circular 
symmetric. 

The subtlety of this theorem is that the autocorrelation is not of statistical circular sym- 
metry only (this is easy to see) but of perfect circular symmetry. The rotated measure R.fi 
is defined by R.fi(A) = n(R~' 1 A) for all measurable sets A. Then fi is of circular symmetry if 
R.fi = [i for all R. 

Proof. Since T is of statistical circular symmetry, T contains tiles Tj congruent to some 
prototile Pi for (countably) infinitely many distinct angles ay = <p(Tj). By equidistribution, 
the frequency of tiles in one certain orientation is 0. But, again by the equidistribution of 
(aj )j>o, the frequency of tiles in a certain range of orientations — say, \J3, 7[C [0, 2tt[ — equals 
the frequency of tiles in orientations \(3 + c, 7 + c[ mod 2tt for all c. The same is obviously 
true for supertiles in T. 

The autocorrelation is defined via difference vectors x— y. We deduce the circular symmetry 
of the autocorrelation by assigning these vectors x — y to supertiles in T . Each such vector 
arises from a constellation of two tiles (usually not adjacent). Consider only vectors x — y 
of fixed a length r > 0. Let k > 1. Then either x and y are contained in the same k-th 
order supertile, or in two different k-th. order supertiles. With growing k, the ratio of the 
latter kind to the former kind tends to zero. Thus it suffices to regard the former kind only. 
Each one arises from a k-th. order supertile, these are equidistributed, hence equidistribution 
of orientations of difference vectors x — y of radius r follows. In other words, the frequency of 
difference vectors x — y of length r with orientation ip £ [/?, 7[C [0, 2ir[ equals the frequency 
of difference vectors of length r in orientations [/3 + c, 7 + c[ mod 2ir for all c. It follows that 
7a = -R-7A f° r all rotations R. □ 

Now we can apply the following standard result, which can be found for instance in 

Lemma 3.2. Let R be an orthogonal map, \i a measure, and let the measure R.fi be given by 
R.H(A) = ^(R^A). Then R.^l = lQl. 

Since i?-7A = 7A, with A as in 13.11 it follows R.^j\ = R.^f\ = ta for all rotations R. 
Altogether we have obtained the following result. 

Theorem 3.3. Let A be as in Theorem \3.1\ arising from a tiling with statistical circular 
symmetry. Then the diffraction of A is circular symmetric. □ 

An immediate consequence is that the pure point part of the diffraction spectrum of a 
tiling with statistical circular symmetry is at most one peak at the origin. A simple argument 
shows that this peak always exists, compare [9], [12]. These results are relevant for physics 
because such tilings can serve as a simple model for powder diffraction [3]. 

4. Tiling dynamical systems 

An important tool in investigating properties of nonperiodic tilings is the fact that such a 
tiling gives rise to a dynamical system, compare [15 ,[20j,[18j,[l], [E]. From the properties of 



6 



D. FRETTLOH 



the dynamical system one can deduce properties of the tiling. For brevity, we describe the 
metrical dynamical system corresponding to [14]. In the context of this paper, the topology 
defined by this metric is the same as the more general local rubber topology in [3], as well as 
the topology in [15J and the local topology in [18]. Let d be the following metric, measuring 
whether two tilings are 'close' to each other. 

d(T, T') = min{-^, inf{e : B 1/e n (T + s) = B 1/e n (R a T' + t), \\s\\, \\t\\ < |, \a\ < e}\ 
y2 e>o 2 

If d(T,T') is small, then T and T agree on a large ball around the origin, after a small 
rotation followed by a small translation. Now, the hull X7- of the tiling T is the closure of 
{T + t : t £ I 2 } with respect to d. It is a well-known fact that the hull of a primitive 
substitution tiling with substitution a is the family X CT of all substitution tilings (compare 
Def. II. ip . Therefore, we denote the hull of T by X<j, if it is clear which substitution a belongs 
to T. 

A substitution a gives rise to a dynamical system (X CT ,G), where either G = M? (all 
translations in the plane), or G = E{2) (all direct Euclidean motions, i.e., all maps x \— ► Rx+t, 
R a rotation, t a translation vector). The study of such dynamical systems yields deep insight 
into the diffraction properties of the tilings in X CT . A milestone in the mathematical diffraction 
theory is the result that the diffraction spectrum is pure point if and only if the dynamical 
spectrum is pure point [9], [20], [18], [11], [3J. Other important results are collected in the 
following theorem. It summarises contributions of several people to the subject over the last 
two decades. For brevity, it is formulated for the special case of plane tilings only. 

Theorem 4.1. Let a be a primitive substitution. 

(1) X(j has uniform cluster frequencies, that is, each patch occurs with a well defined 
frequency throughout the tilings in X CT . Here, frequency is defined analogously to ([!]). 

(2) If the tiles in tilings in X CT show finitely many orientations, then (X CT ,R 2 ) is uniquely 
ergodic. 

(3) If the tiles in tilings in X CT show infinitely many orientations, then (X CT ,G) is uniquely 
ergodic, where G G {R 2 ,E(2)}. 

(4) All tilings in X CT have the same autocorrelation, thus the same diffraction spectrum. 

The first part is a simple consequence of Def. [TTTl compare the discussion in Section [21 or [20] 
for 'self-affine' tilings. (Each X CT as above contains a self-afnne tiling, and by primitivity, all 
tilings in X CT have the same cluster frequencies.) However, for tilings with statistical circular 
symmetry, these frequencies are zero. But see [6] for a similar uniformity result in this case. 
A detained proof of the second part can be found in |20j or [11]. For tilings with tiles in 
infinitely many orientations, part three is stated in [14J for G = E(2). For G = M 2 , it can 
hopefully be shown via Theorem 12.21 adapting methods from [11]. We aim to give a precise 
and detained proof in the future. The last part of the theorem is an immediate consequence 
from part one and the definition of the autocorrelation. 
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